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We study the long-time dynamics of a forest-fire model with deterministic tree growth and in- 
stantaneous burning of entire forests by stochastic lightning strikes. Asymptotically the system 
organizes into a coarsening self-similar mosaic of synchronized patches within which trees regrow 
and burn simultaneously. We show that the average patch length (L) grows linearly with time as 
t — » oo. The number density of patches of length L, N(L,t), scales as {L)~ 2 Af(L/(L}), and within 
a mean-field rate equation description we find that this scaling function decays as Af(x) ~ e -1 ^ for 
x — > 0, and as e~ x for x — » oo. In one dimension, we develop an event-driven cluster algorithm to 
study the asymptotic behavior of large systems. Our numerical results are consistent with mean-field 
predictions for patch coarsening. 
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I. INTRODUCTION 

Forest-fire models ||, [|, ||, Q are simple archetypical 
examples of driven dissipative systems that exhibit in- 
triguingly rich spatio-temporal structures j| [|, |l f?l || [jj . 
These models provide a simple paradigm for cooperative 
time-dependent phenomena, such as epidemics, oscilla- 
tory chemical reactions, electrical neuron activity, car- 
diac dynamics, and turbulence ^ O, [l^, [l3|]. As 
the name indicates, the forest-fire model ostensibly de- 
scribes the time evolution of burning trees in a forest. In 
typical models of this genre, trees are located at regu- 
lar lattice sites and each can exist in one of three states: 
burnt, alive, or burning. The dynamics involves the fol- 
lowing elements: (i) A burnt tree turns into a living tree 
at some specified rate; (ii) a living tree can be ignited, 
either by a lightning strike or by fire spreading from a 
neighboring burning tree; (iii) after a specified time in- 
terval a burning tree is consumed and the fire at this 
location is extinguished. 

Depending on which of these processes are operative, 
as well as their relative rates, different dynamical be- 
haviors can arise, ranging from self-organized critical be- 
havior with fires of all sizes occurring j|, ||, ||], to spi- 
ral fire- front wave propagation [l4j. While forest fire 
models have been extensively investigated, there is still 
uncertainty about their long-time properties even after 
extensive numerical simulations in many realizations of 
the model |§. 

This work is focused on a specific version of the forest- 
fire model which exhibits coarsening [ fl5| rather than self- 
organized criticality or complex fire front propagation. 
Because of this phenomenological simplicity, we can ap- 
ply the rate equations in a natural way to determine the 
evolution of the system. The model itself was first in- 
troduced by Drossel flfj. Its crucial feature is that tree 
growth is deterministic; a tree that has just burned re- 
mains dead for exactly one time unit and then a new tree 
reappears. This particular regrowth rule is the mecha- 
nism that gives rise to a coarsening mosaic of growing 



synchronized forests. This is a generic feature and does 
not require the tuning of model parameters to critical 
values. 

Let us define a "patch" as a coherent region of the 
system that is either occupied by live trees or by burnt 
trees. This patch evolves by tree regrowth and by the 
burning of trees due to lightning strikes. The mecha- 
nism for coarsening is that adjacent patches must even- 
tually synchronize, after which they evolve in phase p6| . 
While neighboring patches begin their existence as dis- 
tinct, eventually the burnt trees in one patch will regen- 
erate while the adjacent patch is still forested. When this 
occurs, all the trees in these two patches become incor- 
porated into an augmented patch which then evolves as 
a single unit. (Fig. [j]). 

In Ref. the process was investigated numerically 
and the total number of patches N(t) was found to de- 
cay with time. However, more quantitative observations 
were not reported. We will show that N(t) oc t^ 1 and 
we will investigate the patch length distribution, both 
analytically and numerically. For our analytic study, we 
will employ the classical rate equation of aggregation ki- 
netics [|17|. This approach is ideally suited to treat the 
coarsening behavior of the system under investigation. 

Another important feature of our approach is that we 
treat the dynamics at a mesoscopic level in which the 
basic units are synchronized patches rather than indi- 
vidual trees. In the long-time limit, we will argue that 
the lifetime of such patches can be viewed as determin- 
istic. By tracking only the merging events of adjacent 
patches we are able to investigate systems of effectively 
much larger size and to much longer times than that ac- 
cessible by tree-based simulations. In addition to obtain- 
ing a power-law growth of the average patch length, our 
method yields clean results for the probability distribu- 
tion of patch lengths. This distribution is found to obey 
scaling, with no memory of the initial state retained and 
with the asymptotes of the distribution in good agree- 
ment with rate equation predictions. While we focus on 
the particular case of one dimension, our approach should 
also apply in higher dimensions. 



In the next section, we define the model and outline the 
effective mesoscopic picture for the evolution of patches. 
In Sec. Ill, a rate equation description for this evolution 
is presented and basic results about the patch length dis- 
tribution are obtained. We explain our simulational ap- 
proach and describe the results that follow in Sec. IV. 
Our basic conclusions are given in the last section. 



II. FOREST EVOLUTION IN ONE DIMENSION 

The evolution of the system is governed by the com- 
petition between two fundamental time scales. Suppose 
that each tree in the system may be struck by lightning 
at a rate A. Then a forest of length L has a characteris- 
tic lifetime t occ ~ (AL) _1 before one of its trees is struck 
by lightning. We assume that the time to burn the for- 
est completely is much less than any other time scale in 
the problem, so that we view the burning of a forest as 
instantaneous. There is also the deterministic time in- 
terval i emp between the instant that a forest burns down 
and the reappearance of trees in this burnt patch. We 
assume that this refractory time is the same for all trees, 
so that regrowth of trees in a single burnt patch occurs 
simultaneously. Without loss of generality, we take this 
refractory time to be t e mp = 1- 

At early times, where t e mp < t occ , the stochastic nature 
of lightning events is important. However, as we shall 
soon show, patches naturally grow with time. Thus t emp 
eventually becomes much larger than t occ and the fire 
dynamics becomes nearly deterministic in the long-time 
limit. To describe this late-stage dynamics, we ignore 
individual trees and treat the system mesoscopically as a 
contiguous array of patches, each with length Lj. Each 
patch can either be forested or burnt. If two distinct but 
adjoining forested patches arise by the regrowth of one 
patch next to a forested patch, they immediately join to 
form a larger patch (Fig. Q). 

At long times, patches have a small lifetime and they 
are almost always in the burnt state. Without loss of 
generality we initialize the system so that it is effectively 
in this long-time state. That is, at t — 0, we consider all 
patches to be burnt, and we define Tj (with < Tj < 1) to 
be the time at which the j th burnt patch first becomes a 
forest. Consider now two adjacent patches, and let <j>j — 
Tj+\—Tj be the time difference between the appearance of 
forest j and forest j + 1. After these two forests undergo 
one cycle of regrowth and subsequent burning, <pj changes 
by Atj+i — Atj, where Atj = (AL_,) _1 is the lifetime of 
the j th forest. This shift in the difference of burning 
times continues until cj>j reaches either <fij = or <pj = 1. 
When this occurs, the forests necessarily join and are 
subsequently synchronized (Fig. [I]). 

Since lightning strikes a forest of length Lj at rate ALj , 
the lifetime Atj of each forest is a stochastic variable 
whose average value is (Atj) = (ALj)' 1 . In the long time 
limit, these lifetimes become very small. Therefore we 
may replace the sum of a large number of these lifetimes 
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FIG. 1: Illustration of the merging of two adjacent patches 
of lengths L\ and L2 > Li. White space indicates a burned 
patch while the shaded region indicates a forested patch. The 
times until the first regrowth events are ri and Ti- Each 
forested patch survives for a time Atj oc LJ 1 until lightning 
strikes and instantaneous burning occurs. The patch then 
regrows after exactly one time unit elapses. Because Ati and 
At2 are different, the two patches will eventually synchronize 
at the joining time Ti^- 



by the average lifetime times the number of cycles. Thus 
for two adjacent patches j and j + 1 that satisfy Lj < 



Lj + i and 



Tj < T 



the joining time Tjj+i is 



T. 
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An analogous result holds when Lj > L 
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these joining times between adjacent patches, we con- 
clude that the typical joining time scales with the aver- 
age patch length as T ~ (At)' 1 ~ A(L). Thus in a time 
interval dt ~ T, a typical patch grows by an amount 
dl ~ L. Consequently 4r(L) ~ T _1 (L), and we obtain 



(L) ~ A" 1 *. 



(2) 



In previous forest-fire models that are driven by 
stochastic tree growth and by stochastic lightning strikes 
H ||, the latter rate needs to be very small to ensure 
non-trivial dynamics. There is no need for such param- 
eter tuning in the present model, as the magnitude of A 
determines only the overall scaling of the typical forest 
length. Therefore, we set A = 1 henceforth. 



III. RATE EQUATION DESCRIPTION 

A natural approach to determine the evolution of the 
patch length distribution is the rate equations. Let 
N(L,t) be the number density of patches of length L 
at time t, and let N(t) = / °° N(L, t) dL be the total 
number of patches of any length. In the rate equation 
description of these quantities, we make the mean-field 
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assumption that there are no correlations between adja- 
cent patches. In a similar spirit, we also ignore the initial 
phase difference between two patches in Eqs. (|l]), so that 
the joining rate of two patches is simply proportional to 
(At,-) — (A< J+ i) oc ILJ 1 — LJ +1 \. With these approxima- 
tions, the patch length distribution evolves according to 
the rate equation 



dN{L,t) 

m 
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2iV(t) Jo 
N(L,t) 
N(t) ,„ 



dl K(l, L-l) N(l, t) N(L - I, t) 
dlK(L,l)N(l,t). (3) 
is the rate at which a patch 



Here K(x, y) — \x~ 1 — y~ 
of length x joins with a patch of length y. 

This rate equation is nearly identical in form to the cor- 
responding equation for the kinetics of aggregation (jl7| , 
except for the overall factor of l/N(t). This difference 
arises because we consider a finite system and we track 
the number of forests of a given length rather than the 
corresponding probability. However, this factor can be 
absorbed into a rescaled time variable defined by 

T=[*-*- (4) 

Jo N(t'Y W 

to reduce Eq. (||) to the standard form of the rate equa- 
tion for aggregation. This can then be analyzed by well- 
established methods fL7| , 

While the rate equation with reaction rate K(x, y) = 
\x~ x — has not been solved, basic features about 
the long-time solution can be inferred from a scaling ap- 
proach. In general, the long-time behavior of the rate 
equations with homogeneous reaction rates that satisfy: 
(i) K(ax,ay) = a x K(x,y), and (ii) K(x,y) ~ x^y v for 
x <C y (A = (i + v), have been generically classified ac- 
cording to whether \i > 0, // < 0, or /i = |l7], In 
all cases, the asymptotic length distribution approaches 
a scaling form 



N(L,t)^(L(t))-*Af(L/(L(t))), 



(5) 



in which the average patch length grows algebraically 
with time, (L) ~ T 1 /( 1_A \ when A < 1. However, 
the scaling function J\f exhibits different behaviors in the 
three cases. 

The reaction rate for our forest-fire model, K(x,y) — 
\x~ l — y~ l \, is homogeneous with homogeneity exponent 
A = — 1, while /i = — 1 and thus v = 0. Therefore, the 
scaling theory prediction for the average patch length 
becomes (L) ~ T 1 / 2 . Consequently, N{T) ~ (L) _1 ~ 
T -1 / 2 , and from Eq. (Q), we recover (L) ~ t, in agree- 
ment with the qualitative argument preceding Eq. (||). 
According to the general classification scheme of van 
Dongen and Ernst [jl7| [l^], since /j, = — 1, the scaling 
function N(x) should vanish exponentially in the limits 
of small and large x: 



N(x) 



-l/x 



For short patches, this result therefore predicts 
N(L,t) - e- l l L . This behavior can also be established 
directly from Eq. (Q). When L <C t, the gain term 
in the rate equation can generally be ignored. Addi- 
tionally, in this limit the reaction rate K(L, I) simpli- 
fies to L~ x . Hence the density of short patches satisfies 
%£f — —L~ 1 N, which indeed implies the above exponen- 
tial decay. 



IV. SIMULATION RESULTS 

In our simulations, we start the system with a random 
array of patches of lengths {Lj(t = 0)}. As discussed in 
Sec. II, it is asymptotically exact to replace the stochastic 
forest lifetime Atj by its average value (At,) = (Lj)~ l , 
The time for regrowth of a tree is always equal to one. 
Thus the dynamical steps become deterministic and ran- 
domness enters only through the initial conditions. 

For convenience, we assume that each patch is initially 
in the burned state. From these regrowth and burning 
processes, a simulation of forest fires should be based on 
the following steps: 

1. Initialize the line with patches of random lengths 
Lj and with random times Tj, for j — 1, 2, . . . , N at 
which the j th patch first turns into a forest. Assign 
a lifetime (Lj)^ 1 to a forest of length Lj. 



2. Use Eq. (|lj), compute the joining times Tjj + i 
all neighboring pairs of patches. 



for 



x 
x 



0, 

00. 



(6) 



3. Sort the list of joining times {Tjj+i} in ascending 
order. A standard sort algorithm |19| requires a 
time of the order of TV In A for a set of N elements. 

4. Join the pair of patches (j, j + 1) with the minimal 
joining time and increment the time accordingly. 
Recompute the joining times Tj-ij and Tj t j+2 of 
the patches adjacent to the newly-joined forest. 

5. Decrement the total number of patches by 1 and 
return to step 3. 

Such an algorithm is perforce inefficient because of the 
re-sorting of joining times after each event. However, 
this step is typically unnecessary for two reasons. First 
we "cache" only a small subset of the joining times with 
Tjj+i less than a judiciously-chosen cutoff time T c and 
sort only this subset in step 3. We need not consider 
joining times 2^ 3+1 > T c because these events in the 
far future will never be considered before the current list 
of joining times needs to be re-sorted. This restriction 
significantly reduces the time needed to sort the joining 
time list for large N. 

Second, it is unnecessary to re-sort this reduced join- 
ing time list after each joining event because only the two 
joining times Tj-i.j and Tjj + 2 are modified. If these up- 
dated joining times are greater than the elements in the 
already-sorted list, there is no need for re-sorting. One 



4 



can simply continue to use the joining times from the 
pre-sorted list to define joining events until one of the 
newly-created joining times becomes less than the next 
joining time in the pre-sorted list. Only when such a 
mis-ordering occurs is it necessary to return to step 3 
and re-sort the joining time list. These steps are com- 
pletely analogous to those employed in Ref. |^0| to simu- 
late the kinetics of one-dimensional ballistic annihilation 
reactions efficiently. 
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FIG. 2: Average patch length (L(t)) versus time t for ini- 
tial length distribution N(L, t = 0) = L^ 1 exp (— L/Lq), with 
Lq — 0.1. A straight line of slope 1 is shown for comparison. 
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FIG. 3: Dependence of the normalized correlation function 
{{LiLi + k) — (L) 2 )/<7 2 , where a 2 — {L 2 } — (L) 2 , on patch sep- 
aration k at t = 2.25 9 w 1478. 

We initialize the system with N = 5 x 10 6 patches 
whose lengths are randomly drawn from the distribution 
N(L,t = 0) = i " 1 cxp(- J L/L ), with L = 0.1. This 
initial length should be viewed as much larger than the 
lattice spacing between individual trees. We are inter- 
ested in the intermediate asymptotic regime, where the 
average patch length is growing systematically with time 
and before finite-size effects begin to play a role. Fig. |2| 




FIG. 4: Number of patches N(L, t) of length L at t = 1478, 
3325, 7482, and 16384 (t = 2.25™ with n = 9 - 12) plotted 
in scaled form. The main plot shows the data on a semi- 
logarithmic scale to illustrate the exponential decay of the 
large-length tail. The inset shows the same data in the small- 
length limit on a linear scale to highlight the essentially sin- 
gular form. 



indicates that this intermediate asymptotic regime be- 
gins when t > t* sw 10. An important feature of the sys- 
tem at long times is that there is only a very short-range 
spatial correlation in the lengths of neighboring patches 
(Fig. ||). In particular, the normalized correlation func- 
tion {(LiL i+k ) - (L) 2 )/a 2 , with a 2 = (L 2 ) - (L) 2 , quickly 
approaches zero for k > 2. This provides empirical justi- 
fication for the validity of the mean-field approximation 
of the rate equations. It is worth remarking that such 
lack of spatial correlations appear ubiquitously in many 
one-dimensional coarsening processes [fl5| [jlf . 

We now examine the behavior of N(L,t) for represen- 
tative values of t > t* (Fig. ||) . The different sets exhibit 
data collapse according to the scaling ansatz of Eq. (0). 
The large-length tail of this distribution appears to be 
consistent with a simple exponential decay. However, 
there is a very small but apparently systematic downward 
curvature to the data for which we do not have an expla- 
nation. The small-length tail decays extremely rapidly 
near the origin and there are essentially no patches with 
scaled length less than 0.2. This is consistent with the 
essential singularity predicted by Eq. (|^). Again, it is 
worth remarking that in almost all coarsening processes 
that are controlled by an underlying diffusion process, the 
small-length tail of the patch-length distribution is a lin- 
ear function. The one well-known example of an essential 
singularity in the small-size tail of a cluster distribution 
is the aggregation of Brownian particles M| . 



V. CONCLUSIONS 

We have developed a mesoscopic description for a 
forest-fire model with stochastic lightning strikes and de- 
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terministic tree growth. Instead of treating the system at 
the level of single trees, the basic element in our descrip- 
tion are patches of synchronized trees. Each patch un- 
dergoes periodic cycles of burning and regrowth, and in 
the long-time limit, the lifetime of patches can be viewed 
as deterministic. Whenever two adjoining patches are si- 
multaneously in the forested state, the patches join and 
remain synchronized forever. The number of distinct 
patches decreases while their typical length grows con- 
tinuously with time as in classical coarsening processes 

This mesoscopic description is well-suited to a rate 
equation approach for the evolution of the patch-length 
distribution, as well as efficient simulations. Numerically, 
we find that the average patch length grows linearly with 
time, while the number of patches correspondingly de- 
creases as 1/t. The patch length distribution is sharply 
peaked, with an exponential large-length tail and an es- 
sential singularity in the small-length tail. These fea- 
tures are consistent with the rate equation predictions. 
Even though the system has nearest-neighbor interac- 
tions only, there are essentially no spatial correlations in 
the lengths of neighboring patches. Because of this lack of 
spatial correlation, we can anticipate that the rate equa- 
tion predictions, which are based on no correlations be- 
tween patch lengths, should provide an accurate account 
of the one-dimensional simulations. 

Most aspects of our approach can be extended to 
higher spatial dimensions d. A complicating factor in 



developing numerical simulations is that the number of 
neighbors for a given patch is variable. Nevertheless, the 
same updating rule given by Eqs. (Q) will still apply, with 
patch length being replaced by patch volume. As a result, 
we expect that the average patch volume should grow 
linearly with time. Under the assumption that patches 
remain compact, this would imply that the typical length 
scale of a growing patch would grow in time as L ~ t x l d . 
It is of course far from obvious that patches remain com- 
pact. Understanding of the forest-fire model in higher 
dimensions appears to be an interesting challenge. 

Finally, we want to stress that the seemingly innocent 
change of the tree growth rule from stochastic to deter- 
ministic drastically affects the dynamics. Almost all ear- 
lier work focused on stochastic tree growth. According 
to the literature on such model, self-organized critical be- 
havior should occur in the limit of infinitesimal rate of 
lightning. However, after considerable effort, the under- 
standing of the scaling laws is still incomplete, and even 
their very existence has recently been questioned |9|. In 
contrast, when the tree growth is deterministic, there is 
no need to tune the rate of lightning strikes to zero, and 
the model exhibits simple coarsening rather than com- 
plex time-dependent phenomena. It is of course possible 
that the higher-dimensional version of the model will of- 
fer some surprises. 

We are grateful to NSF grant DMR9978902 for partial 
support of this work. 
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